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The scattering of particles by a Coulomb field is an important subject in quantum mechanics 
as well as in classical mechanics, because the Coulomb field describes real physical interations 
between charged particles, say, electrons and nucleuses. In the nonrelativistic quantum me- 
chanics, the problem is described by a Schrodinger equation. In ordinary three-dimensional 
space, this equation can be solved in rotational parabolic coordinates, and the scattering am- 
plitude is obtained in a closed form. As a special case of a centrally symmetric field, the 
Schrodinger equation can also be solved in the spherical coordinates, and the scattering am- 
plitude is obtained as a partial wave series. This has been well known for a long time, and 
the solutions can be found in the popular textbooks of quantum mechanics [1,2]. However, it 
appears that the partial wave amplitude is not evaluated in the textbooks. Were it not for the 
solution in the parabolic coordinates, how can one directly sum up the partial wave series to 
obtain the closed form? A solution can be found in an early paper of Mott [3] . However, the 
calculation is not well known since it is not given even in the famous monograph on collision 
theory by that author and Massey [4]. Other books on scattering theory do not deal with the 
matter either [5,6]. The purpose of this note is to evaluate the partial wave series in a simpler 
way. The method is easy for the students and thus may be of some interest for the pedagogical 
purpose. 

Let us write the Hamiltonian in the form 



where fj, is the reduced mass of the system, and k > (k < 0) corresponds to an attractive 
(repulsive) force. The partial wave scattering amplitude obtained by solving the problem in 
spherical coordinates reads [1,2] 
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where k = y/2 fiE/h (E > is the energy of the incident particle), and 
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where 5i is the phase shift of the /th partial wave, and (3 = fiK/h 2 k = n/hv where v is the 
incident velocity. Since [1,2] 

oo 

£(2Z + l)P,(cos0) = 4<S(cos0-l), (4) 

1=0 

we have for 6 ^ 

1 oo 

/(0) = ^E(2J + l)Wcos0). (5) 

We define 

oo 

</(x)=X;(2/ + l)5,P,(x), (6) 
z=o 

and are going to evaluate it. Using the relation 

(2l + l)P l (x) = P; +1 (x)-PUx) } (7) 
where P-i(x) = 0, and define an auxiliary function 

oo 

G(x) = J2S l [P l+1 (x)-P l ^ 1 (x)}, (8) 

we have 

g(x)=G'(x). (9) 

Next using the ralation 

(2/ + l^Ptix) = (/ + l)P m (x) + LFU(:r), (10) 
and the definition (8), we have 

oo 

xg(x) = £S,[(Z + 1 - + + i/3)Pi-i(x)] + if3G(x). (11) 

1=0 



Then making use of the relation 



T(z + l) = zT(z), (12) 
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and Eq. (3), we have 

(Z + 1 - 0)Si = (I + 1 + (Z + i0)Si = (Z - (13) 

Substituting these relations into Eq. (11), after some simple algebras we have 

xg(x) = g(x) + i(3G(x) - i(3S . (14) 

Taking Eq. (9) into account, we obtain 

(1 - x)G\x) + i(3G{x) = if3S . (15) 

This is an ordinary differential equation for the auxiliary function G(x). It can be easily solved 
to give 

G(x)=Aexp[ip\n(l-x)]+S , (16) 

where A is a constant of integration. To determine A we consider the series (8) at x — — 1. 
Note that Pi(— 1) = (—1)' and = 0, only the term with I = contributes and so 

G(—l) = —S . Substituting this into the above result one can find A and arrives at 

G(x) = -2S exp[i/3m((l-x)/2)] + S . (17) 

It should be remarked that in getting Eq. (9) we have differentiated G(x) term by term. This 
is not allowed at x — 1 since g(x) is obviously singular at the point. Thus Eq. (15) does not 
hold at x — 1, and it is not surprising to find that the result (17) does not agree with the 
definition (8) at that point. By using Eq. (9) we can easily find g(x). Setting x = cos9 and 
multiplying it by the coefficient l/2ik we arrive at 

r(l-^)exp^lnsin 2 (g/2)] 
I{} iV{i(3) 2ksm 2 (6/2) ' 1 j 

This is the closed result obtained in the parabolic coordinates [1,2]. Two remarks about the 
result. First, we have dropped a 5(cos9 — 1) term for 9 ^ in the above calculation. Here the 



final result is more singular than S(cos9 — 1) when 9 — > 0, thus that term can indeed be dropped 
everywhere. Second, it has been argued that f{9) in the form of Eq. (2) or (5) is divergent 
even when 9 ^ [7]. The final result (18) is well defined everywhere except at 9 — 0, however. 
Therefore the above calculation may also be regarded as a regularization procedure. (The 
main point lies in the term by term differentiation pointed out above, which is mathematically 
not allowed even when i / 1.) The validity of this procedure is obvious because the final 
result (18) coincides with the one obtained from exact solutions in the parabolic coordinates. 
It should be pointed out that the calculation of Mott [3], though rather different from ours, 
involves some similar story. 

Mathematically one can consider the Coulomb scattering problem in a general space di- 
mension d > 2. In d = 2 the problem has been solved in both parabolic coordinates [8-9] and 
polar coordinates [10]. The scattering amplitude is obtained in a closed form in the former 
case and as a partial wave series in the latter case. The method for summing up the partial 
wave series in d = 2 is rather different from that in d = 3 [10,11]. In more general dimensions 
d > 3 the problem has been solved in parabolic coordinates and a closed form for the scattering 
amplitude was obtained [12]. If one solves the problem in d- dimensional spherical coordinates, 
one would obtain for the scattering amplitude a partial wave series similar to Eq. (2) , where 
the Legendre polynomials are replaced by the Gegenbauer ones. We expect that the method 
given by Mott or the one given here may be efficient for summing up that series as well. 
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munications and encouragement. This work was supported by the National Natural Science 
Foundation of China. 
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